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1. Introduction 

An important dynamical mechanism, by which fundamental particles acquire masses, is 
through interactions with vacuum condensates. Thus, a massive particle may carve out 
around itself a spherical region P] or a shell [2j in which the condensate is suppressed, 
thus reducing the effective mass of the particle at the expense of volume and gradient 
energy associated with the condensate. This picture has interesting phenomenological 
consequences jU 0] . 

This dynamical distortion of the homogeneous vacuum condensate configuration, 
namely, formation of fermion bag solitons, was demonstrated explicitly by Dashen, 
Hasslacher and Neveu (DHN) |lj many years ago, who studied semiclassical bound 
states in the 1-1-1 dimensional Gross-Neveu (GN) model using the inverse scattering 
method jHI- Following DHN, Shei [Zj has applied the inverse scattering method to study 
solitons in the 1-1-1 dimensional Nambu-Jona-Lasinio (NJL) model jH] in the large-iV 
limit. 

Fermion bags in the GN model were discussed in the literature several other times 
since the work of DHN, using alternative methods [HlCniE]- For a review on these and 
related matters (with an emphasis on the relativistic Hartree-Fock approximation) see 
. For a more recent review of static fermion bags in the GN model (with an emphasis 
on reflectionless backgrounds and supersymmetric quantum mechanics) see jlH]. The 
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large- semiclassical DHN spectrum of these fermion bags turns out to be essentially 
correct also for finite A^, as analysis of the exact factorizable S-matrix of the GN model 
reveals [T^ . 

A variational calculation of these effects in the 1 + 1 dimensional massive 
generalization of the Gross- Neveu model, which we will refer to as MGN, was carried in 
[T3] a few years ago, and more recently in [16^,. Very recently, we studied static fermion 
bags in the MGN model ^Tj, which we obtained using the inverse-scattering formalism, 
thus avoiding the need to choose a trial variational field configuration. Our main result 
in jTTj is that the extremal static bag configurations are refiectionless, as in the massless 
Gross- Neveu model. In the next section we briefly review the results of jEj, leaving 
technical details out. Then, in Section 3, we show that a subclass of the reflectionless 
solitons of arise self-consistently in the 1 + 1 dimensional massive NJL (MNJL) 
model. The latter extends the results of (3 ^] for the massless NJL model. Solitons in 
the MNJL model were also recently studied in ^HI, where a derivative expansion was 
carried out around a particular soliton background of the corresponding massless NJL 
model. 

2. Solitons in the Massive Gross-Neveu Model 

One way of writing the action for the MGN model is 

^ = j d^^\^^i^a[i0-<y]^a-i^A<y^ -^Ma)^ , (1) 

where the ipa {a = 1, . . . , A^,) are massive Dirac fermions and a is an auxiliary fleld. 
Integrating the a out results in an equivalent form of (^, with quartic fermion self- 
interactions. 

An obvious symmetry of (0) with its A^ Dirac spinors is U{N). Actually, is 
symmetric under the larger group 0{2N) jlj (see also Section 1 of [iH])- The fact that 
the symmetry group of ^ is 0{2N) rather than U{N) is related to the fact that it is 
invariant against charge- conjugation, like the massless GN model. Consequently, the 
energy eigenvalues of the Dirac equation associated with (jT}, [i^ — cr{x)] 4> = 0, come in 
±uj pairs. 

As usual, the theory (0) can be rewritten with the help of the scalar flavor singlet 
auxiliary fleld cr(x). Also as usual, we take the large A^ limit holding A = Ng"^ flxed. 
Integrating out the fermions, we obtain the bare effective action 

S[a] = — ^ [ d^x [a^ - 2Ma) - iN Tr log {i$ - a) . (2) 

Noting that '^^{i$ — a) = —{i0 + c'")75, we can rewrite the Tr log{i0 — a) as 
|Tr log — a){i0 + cr)] . In this paper we focus on static soliton conflgurations. If a 

is time independent, the latter expression may be further simplifled to |^ J i^lTr log(/i+ — 
tu^) + Tr \og{h_ — uj"^)] where h± = —d"^ + a"^ ± a' , and where T is an infra-red 
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temporal regulator. As it turns out, the two Schrodinger operators h± are isospectral 
(see Appendix A of [13J and Section 2 of ^1]) and thus we obtain 



oo 



S[a] = j d^x {a^ -2Ma) -tNT j ^Ti \og{h_ - u^) . (3) 

— oo 

In contrast to the standard massless GN model, the MGN model studied here is not 
invariant under the Z2 symmetry ip Ts^/', cr — > —a, and the physics is correspondingly 
quite different. As a result of the Z2 degeneracy of its vacuum, the GN model contains 
a soliton (the so called CCGZ kink 01111311111131) in which the a field takes on equal 
and opposite values at x = ±00. The stability of this soliton is obviously guaranteed by 
topological considerations. With any non-zero M the vacuum value of a is unique and 
the CCGZ kink becomes infinitely massive and disappears. If any soliton exists at all, 
its stability has to depend on the energetics of trapping fermions. 

Let us briefly recall the computation of the unique vacuum of the MGN model. 
We shall follow ^H]- For an earlier analysis of the MGN ground state (as well as its 
thermodynamics), see [201 • Setting cr to a constant we obtain from Q the renormalized 

2 2 1 2 

effective potential (per flavor) V{a,fi) = ^ log ^ + ^ — M{^)a , where /i is a 
sliding renormalization scale with A(/i) = Ng'^{fi) and M{fi) the running couplings. By 
equating the coefficient of o"^ in two versions of V, one defined with fii and the other 
with fi2, we find immediately that — = ^ log ^ and thus the coupling A is 
asymptotically free, just as in the GN model. Furthermore, by equating the coefficient 
of cr in y we see that the ratio is a renormalization group invariant. Thus, M and 
A have the same scale dependence. 

Without loss of generality we assume that M(/i) > and thus the absolute 
minimum of V{a,fi), namely, the vacuum condensate m = (a), is the unique (and 
positive) solution of the gap equation ^ = m ^ log ^ + — ^^j^ = . 

Referring to (PJ, we see that m is the mass of the fermion. Using the explicit scale 
dependence of A(/i), we can re- write the gap equation as = which shows 

manifestly that m, an observable physical quantity, is a renormalization group invariant. 
This equation also implies that M(m) = m, which makes sense physically. 

Fermion bags correspond to inhomogeneous solutions of the saddle-point equation 
j^rj^ = 0. In particular, static bags cr(x) are the extremal configurations of the energy 
functional (per flavor) £^[cr(x)] = —^^jj^ , subjected to the boundary condition that 
a{x) relaxes to its unique vacuum expectation value m at x = ±00. More specifically, 
we have to evaluate the energy functional of a static configuration cr(x), obeying the 
appropriate boundary conditions at spatial infinity, which supports K pairs of bound 
states of the Dirac equation at energies ±a;„, n = 1, . . . ,K (where, of course, u"^ < m?). 
The bound states at ±ci;„ are to be considered together, due to the charge conjugation 
invariance of the GN model. Due to Pauli's exclusion principle, we can populate each 
of the bound states itUn with up to fermions. In such a typical multiparticle state, 
the negative frequency state is populated hj N — hn fermions (i.e., by hn holes or 
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antiparticles) and the positive frequency state contains Pn fermions (or particles). We 
shall refer to the total number of particles and antiparticles trapped in the n-th pair of 
bound states i'n= Pn + hn as the valence, or occupation number of that pair. 

The energy functional £^[(T(a;)] is, in principle, a complicated and generally unknown 
functional of (t{x) and of its derivatives (which furthermore, requires regularization). 
Thus, the extremum condition = 0, as a functional equation for a{x), seems 

intractable. The considerable complexity of the functional equations that determine the 
extremal cr(a;) configurations is the source of all difficulties that arise in any attempt to 
solve the model under consideration. DHN found a way around this difficulty in the case 
of the GN model They have used inverse scattering techniques jH] to express the 
(regulated) energy functional S[a] in terms of the so-called "scattering data" associated 
with, e.g., the hamiltonian mentioned above (and thus with cr(x)), and then solved 
the extremum condition |^|^ = with respect to those data. 

The scattering data associated with /i_ are [H] the reflection amplitude r{k) of the 
Schrodinger operator /i_ at momentum k, the number K of bound states in h- and their 
corresponding energies < tu^ < , {n = 1, . . . , K), and also additional K parameters 
{c„}, where c„ has to do with the normalization of the nth bound state wave function ipn 
of h_. More precisely, the nth bound state wave function, with energy u"^, must decay 
as ipn{x) ~ const, where < = — cu"^ . If we impose that 

ipn{x) be normalized, this will determine the constant coefficient as c„. (With no loss 
of generality, we may take c„ > 0.) Recall that r{—k) = r*{k), since the Schrodinger 
potential V{x) = cr^(x) — (t'{x) is real. Thus, only the values of r{k) for k > enter 
the scattering data. The scattering data are independent variables, which determine 
V{x) uniquely, assuming V{x) belongs to a certain class of potentials which fall-off fast 
enough toward infinity. (Since the MGN does not bear topological solitons, neither h_ 
nor h+ can have a normalizable zero energy eigenstate. Thus, all the Un are strictly 
positive.) 

We can apply directly the results of DHN in order to write down that part of £^[cr(x)] 
which is common to the MGN and GN models, i.e., £^[(T(a;)] with its term proportional to 
M removed, in terms of the scattering data. For lack of space we shall not write DHN's 
expression for the energy functional explicitly. Suffice it is to mention at this point that 
the "DHN-part" of 6^[(T(a:)] depends on the reflection amplitude only via certain regular 
dispersion integrals of the quantity log[l — |r(A;)p]. The well-known reflectionless nature 
of solitons in the GN model is a direct consequence of this simple fact. 

In order to complete the task of expressing the effective action of the MGN model in 
terms of the scattering data, we have to flnd such a representation for the remaining piece 
of £^[(t(x)] proportional to M, namely, for — ^ (c"(2^) "~ '^) dx. The latter integral 
cannot be expressed in terms of the scattering data based on the trace identities of 
the Schrodinger operator /i_ discussed in Appendix B of jl]. Evidently, new analysis is 
required to obtain its representation in terms of the scattering data. Happily enough. 
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we were able to obtain such a representation in ^7j, which reads 

1 [ log [1 — |r(A;)p] f m — K 



^ ^ ^ ^ 2m J im-k ^ \m + KnJ 

— oo 

Thus, the M-dependent part of £^[cr(x)], like its "DHN-part", depends on the reflection 
amplitude only via the combination log[l — |r(A;)p]. Combining these two terms together, 
it follows that d£''^3^a{x)]/6r{k) = F{k)r*{k)/{1 - \r{k)\'^), where F{k) is a calculable 
function, which does not vanish identically. Thus, r{k) = is the unique solution of the 
variational equation 6S^'^^[a{x)]/6r{k) = 0. Static extremal bags in the MGN model are 
reflectionless, as their counterparts in the GN model. 

Explicit formulas for reflectionless cr{x) configurations with an arbitrary number 
K of pairs of bound states are displayed in Appendix B of ^3] • In particular, the one 
which supports a single pair of bound states at energies icjf, {k = ^Jm? — lu^), the one 
originally discovered by DHN, is 

/ 1 m + k\ ( 1 m + K 

tanh \KX log — tanh H — log 



o\x) = m+K , , 

4 m — K J \ 4 m — K 

We see that the formidable problem of finding the extremal (t{x) configurations of 
the energy functional S [a] is reduced to the simpler problem of extremizing an ordinary 
function £^(ti;„,c„) = £ [a{x;u!n,Cn)] with respect to the 2K parameters {cn,ujn} that 
determine the reflectionless background a{x). If we solve this ordinary extremum 
problem, we will be able to calculate the mass of the fermion bag. This we did in 
detail in jT2j. Let us sketch the procedure and state the final result: 

The bare regulated energy function £{ujn) which depends on the bare couplings A 
and M and on the UV-cutoff A explicitly can be renormalized, in a manner essentially 
similar to the effective potential, as was described above. £ is independent of the c„'s, 
which appear in the scattering data. (The latter are thus fiat directions for the energy 
function and determine the collective coordinates of the soliton.) The renormalized 
energy function thus obtained is a sum of the form J2n=i fi^n, ^n) where /(a;, v) is a 
known function, which depends also on the physical mass m explicitly, and also through 
the RG-invariant ratio 7 = = -^^^rrr- Thus, the extremum condition fixes each Un 
in terms of the number of the total number of particles and holes trapped in the 
bound states of the Dirac equation at ±0;^, and not by the numbers of trapped particles 
and holes separately (see ©). This fact is a manifestation of the underlying 0{2N) 
symmetry, which treats particles and holes symmetrically. Moreover, it indicates ^7] 
that this pair of bound states gives rise to an 0{2N) antisymmetric tensor multiplet of 
rank i/„ of soliton states. (As it turns out, only tensors of ranks < < correspond 
to viable solitons |IZI-) The soliton as a whole is therefore the tensor product of all these 
antisymmetric tensor multiplets. Finally, we showed in ^7] that only the irreducible 
{K = 1) soliton was protected by energy conservation and 0{2N) symmetry against 
decaying into lighter solitons (or free massive fermions). Its profile is given by (jSJ, 
where k = msinO (or, equivalently, u = mcos6), with < < 7r/2, and where 6 is 
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determined by the extremum condition 

^+T,tau» = ^. (6) 

The left-hand side of dH} is a monotonically increasing function. Therefore, (jHl) has a 
unique solution in the interval [0, 7r/2]. This solution is evidently smaller than 6'^^ = 
the corresponding value of 6 in the GN model for the same occupation number. Thus, 
the corresponding bound state energy uo = m cos 9 in the MGN model is higher than its 
GN counterpart, and thus less bound. The soliton mass (i.e., the renormalized energy 
function, evaluated at the solution of (jHI)) is 

A1(z/) = A^m -sm 6^ + 7 log . (7) 

\ vr 1 — sin 6' / 

This coincide with the corresponding results of variational calculations presented in 

|15| I16j. which were based on (jS)) as a trial configuration. In fact, it was realized in 

that the trial configuration © is an exact solution of the extremum condition |^|^ = 0, 

provided (0) is used to fix /t = m sin 6*. 



3. Refiectionless solitons in the Massive Nambu-Jona-Lasinio Model 

It is natural to inquire whether the results of the previous section carry over to the phe- 
nomenologically interesting MNJL model. The action for the MNJL model may be writ- 
ten as a generalization of (H)), S" = / d?x j^^^ — (cr + «vr75)] "^a — ^ (^"^ + tt^ — 
where 7r(x) is a pseudo-scalar auxiliary field. (Here we assumed that the 2x2 chiral 
mass matrix does not have a pseudo-scalar component, but this does not restrict the 
generality of our discussion in any way. This particular orientation of the mass matrix 
can be always reached at by performing a global - and therefore, anomaly free- chiral 
rotation in the cr — vr plane.) 

As in our discussion of the MGN model, we can integrate out the fermions, and 
obtain the bare effective action 

S[a] = j d^x {a^ + vr^ - 2Ma) - iN Tr log {i0 - a - in-f^) . (8) 

As before, we take the large N limit, holding A = Ng"^ fixed. Unlike the NJL model, 
with its continuum of degenerate vacua, the ground state of the MNJL model (jH} is 
unique, as in the MGN model. It corresponds to a constant field configuration, where 
vr = and where a = m is determined by an equation identical to the one which arises 
in the MGN model. 

Shei [7j has studied static solitons in the NJL model (i.e., M = in (jSJ) using 
inverse scattering techniques. Similarly to DHN's results for the GN model, he concluded 
that extremal soliton profiles are refiectionless. Some of his results were rederived in 
[TH] . using a certain method based on properties of the diagonal resolvent of the Dirac 
operator, (which was applied first to the GN model in [TT]). 

Could Shei's analysis be extended to study solitons in the MNJL model, similarly to 
the extension of DHN's inverse scattering analysis to the MGN model? Could it be that 
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the self-consistent static soliton backgrounds in the MNJL model are reflectionless? 
It seems that all we need in order to answer these questions is a generalization of 
to the case in which the Dirac operator involves a pseudo-scalar background 7r(x). 
Unfortunately, we were not able (so-far) to find such a generalization, and therefore we 
cannot answer these questions in general at the moment. However, we were able to find 
a particular family of self-consistent reflectionless static solitons in the MNJL model by 
making an educated guess, as we shall now explain. 

The spectrum of the Dirac equation associated with (jH} is not invariant against 
charge-conjugation, unless 7r(x) = 0. Thus, the bound states corresponding to a static 
soliton background are not paired, in general. In particular, as has been shown by Shei, 
there exist solitons in the NJL model which bind fermions into a single bound state. 
However, he has also found solitons with charge-conjugation-invariant spectrum (see 
Eqs.(3.25)-(3.28) in [7|), with a pair of bound states io;;,, in which 7r(x) = identically, 
and (j{x) is given by (0), which thus coincide with the DHN solitons in the GN model, 
for which ujb = mcos (^) and A^dhn('^) = ^^^sin (^). However, unlike in the GN 
model, in the NJL model we must choose p = h = ^ (i.e., a soliton of this type must 
trap an equal number of fermions and anti- fermions) . The reason for this restriction is 
not hard to understand physically: the total chiral rotation A6, namely, the difference 
in arctan between the two ends of the one dimensional space, must be related to the 

fermion number charge Uf deposited in the soliton according to A^^ = — |2I] (see 
also Eqs. (5.10) and (5.22) in jTH|). The soliton profile ((t(x), 7r(x)) under consideration, 
starts at the vacuum point (m, 0) at a; = — oo and returns to it at x = +oo. Thus, 
A6 = nf=p — h = for this soliton. 

Now, any static soliton profile in the MNJL model must start at the unique 
vacuum (m, 0) at x = — oo and return to it at x = -|-oo. Thus, it should bring about 
null total chiral rotation, precisely as Shei's charge- conjugation-invariant configuration 
does. Therefore, if the MNJL bears reflectionless static solitons, they must be of this 
form (or charge-conjugate invariant generalizations thereof, with more pairs of paired 
bound states). The only thing that should change compared to the NJL model is the 
quantization condition, relating Ub and u. 

We have verified that this is indeed the case, simply by substituting this 
configuration into the static inhomogeneous saddle-point equations associated with (jHl). 
Varying (jSJ with respect to 7i{x) we obviously obtain an equation identical in form to 
that of the NJL model. (For the latter, see the second equation in (5.1) in [IB])- Using 
the explicit expressions for the entries of the diagonal resolvent of the Dirac operator 
with a reflectionless {a{x),n{x)) background with two bound states (Eqs. (4.13) and 
(4.14) in [IHj with paired bound states uj2 = —uJi), we see that n{x) = is indeed a 
solution of that equation. (Here, having p = /i = | is essential.). This vr-equation 
leaves lji an undetermined function of u. We still have to vary with respect to cr(x). 
Substituting the explicit expressions for the appropriate entries of the diagonal resolvent 
of the Dirac operator (Eqs. (4. 13), (4.14) and (2.10) of [TH]) in the saddle-point equation 
arising from variation with respect to a{x), and using the simplifying identity Eq. (2.24) 
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of [22]; we arrive simply at the static saddle-point equation of the MGN model, which is 
solved by a{x) given by (0) and the quantization condition ®, leading to soliton mass 
(jZj). Thus, a restricted subset of the extremal refiectionless solitons of the MGN model 
appear, not surprisingly, also in the MNJL model. For these solitons 7r(x) = 0. The 
question whether these solitons exhaust all possibilities in the two-dimensional MNJL 
model remains open. 
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